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Abstract 

We give a brief review of recent developments by the variational Monte Carlo method, in 
addition to some new results. We discuss t-J-type models: the ordinary t-J model in one and two 
dimensions, and the one-dimensional supersymmetric t-J model with the long-range coupling of 
inverse square. 
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°3 ■ 1 Introduction 
O 

\Q . The variational Monte Carlo (VMC) method has been extensively developed for many-body systems, 
since the first successful applications to bosons [jl|and fermions. J|] Generally speaking, it was not until 
| the VMC method appeared that the many-body variation theory came to give reliable results without 
uncontrollable approximations. And now the VMC is a standard method to estimate expectation 

S values in this theory. 
In solid state physics, the VMC method was first applied to lattice systems such as the Hubbard 
model ||, [|, H and the periodic Anderson model, || and modified some prior results obtained in the 
additional approximations. Since reviews of the recent VMC studies in this field have already been 
given, 0, |8[ we avoid the overlap of the material, and concentrate on the recent developments for the 



ordinary t-J model || and the long-range t-J model. [ 10 1 



Before going into specific models, we give a brief summary of the VMC method in §2. In §3.1 we 
discuss the one-dimensional (ID) t-J model, for which many exact behaviors are known. Nevertheless, 
we show how further knowledge is added by comparing the VMC results with the exact and other 
results. This is ascribed to the merit of the explicit form of the wave function. In §3.2 we discuss the 
t-J model on the square lattice, comparing with the properties of the ID system. In §4 we develop a 
new usage of the VMC method, namely identification of the eigenstates by checking the variance. We 
apply it to the supersymmetric long-range t-J model with the coupling of inverse square. 

2 Formulation 



In this section, we summarize the fundamentals of the VMC method. [11] This method is a combination 
of the variation theory with the Monte Carlo technique for evaluating expectation values. Following the 
variation theory, one has to estimate the energy expectation value with respect to a given trial function 
■0 7 (R) Jl2| as a function of 7. [13] Here 7 denotes a set of variational parameters and R represents a 
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certain electron configuration. Next, by using the optimal 7, one should calculate physical quantities. 
Actually, the variational expectation value of an operator O is written as, 



{ ° h ~ <^(R)K(R)) } V 7 (R) ' ( j 



where P 7 (R) is the probability density for a configuration R: 



Er'I^(R')I s 



P t(R) = ^ 11 li2 - ( 2 ) 



According to the Monte Carlo method, the summation over R is replaced by the importance 
sampling with weight P 7 (R). Then, eq.(2.1) is reduced to 

in\ 1 °V'7( R m) /on 
(0) --M^^(R^' (3) 

where R m is the m-th sample picked out and M is the total number of samples. Since -P 7 (R) > 0, 
there is no sign problem. In the ordinary usage of the VMC method, one should reduce the statistical 
errors by taking sufficient samples and keeping statistical independence among samples. 

Inversely one can take advantage of this statistical property to identify the eigenstates, j2] as will 
be seen in §4. When one applies the Hamiltonian to ^> 7 (R), the following formula holds: 

H^(R) = £ 7 V 7 (R) + £ 7 (R), (4) 

where £ 7 (R) is orthogonal to ?/> 7 (R) and = (W) 7 . If ?/> 7 (R) is an exact eigenstate, eq.(2.4) has 
to be reduced to the Schrodinger equation, so that £ 7 (R) = 0. This relation and the definition of 
variance, 

.2 _ „ N2 V _ (g 7 (R)|g 7 (R)) 

(^ 7 (R)|V 7 (R)) 



0% = {{h - e,)\ = ;t)^\tz j i (5) 



lead to the fact that there is no statistical fluctuation for an exact eigenstate. [14] Then, E^ becomes 
the eigenenergy. On the other hand, an approximate function results in an appreciable variance in 
proportion to \j\f~M. Thus, in the actual VMC computations, the vanishing of the fluctuation can 
be definitely distinguished. 

Now let us turn to the trial wave function. For t- J-type models two-body Gutzwiller-Jastrow-type 
wave functions, 

* = n n i 1 - u - *?m) n i- n ^\ $ ( 6 ) 

are often used, where r](rje) is a spin-independent correlation factor (r/(0) = 0) with rj£ = \rj — ri\, 
and $ a Hartree-Fock-type wave function. By specifying $ (one-body state) and r] (two-body effect), 
one can describe a variety of states from metallic to ordered or insulating ones by eq.(2.6). 

Leaving individual forms in the following sections, here we discuss only the Gutzwiller wave function 
(GWF) [jl6|] \&G; which is the exact ground state for the long-range t-J model, and is often used as a 
starting trial function for its simplicity. For the GWF, rj(r) = 1 and <!> = which is the Fermi sea. 
The properties of the GWF were investigated for the Hubbard model; 17, 18] the GWF is generally 
metallic, having a discontinuity of momentum distribution n(k) at k = hp. And a repulsive intersite 
correlation factor is needed to represent the Hubbard model properly, p 



3 t-J Model in One and Two Dimensions 

In this section we discuss the t-J model defined as 

H = H t + Hj, (7) 
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Hj = J 53( s « • s i ~ \ n i n i)i 

(ij) 



(8) 



(9) 



in the subspace with no double occupancy of each site with t, J > 0, and (ij) denotes a nearest- 
neighbor pair. Spin operators vanish when they act on empty sites. Henceforth we take t as the unit 
of energy. 

When we regard eq.(3.1) as an effective Hamiltonian of the Hubbard model or the d-p model in 
the strong coupling regime, the following term of the second-order perturbation has to be considered: 



n 3 



' 4 



53 ( C }-cr C j-v C ]+Tcr C j+T'cr + C ]+ T -cr C j-^ C ]cr C J : +r'a ) 



(10) 



where r and r' indicate vectors to the nearest-neighbor sites. When eq.(3.4) is derived from the 
Hubbard model, J3 = J, [jl^, ||] while J3/J < (ferromagnetic) when derived from the d-p model 
with plausible parameters for high-T c superconductors. 



3.1 One- dimensional case 

Many properties in the one-dimensional (ID) model of eq.(3.1) have already been elucidated by 



various exact methods. [|21[ In the following, we would like to show what kind of new knowledge the 
variation theories have added to those. 

First of all, we take the GWF for this system. In Fig.l we show the total energy for the supersym- 
metric case (J/t = 2) as a function of electron density n = N e /N (N e and N being electron and site 



number respectively). [22] Astonishing facts are that the results of the Bethe Ansatz and the GWF 
agree almost perfectly, and the non-interacting system joins for small n. The three results coincide 
in the low density limit up to the order of n 3 , and the difference is about 0.15% for n = 1 between 
the GWF and the exact result. |4|, 15] This good agreement of the two results is not restricted to the 
energy; momentum distribution and correlation functions also show quantitative agreement, except 
for the long-range behaviors 



23] discussed later. 

Intuitively we can understand these results as follows: for the supersymmetric case, the Hamilto- 
nian eq.(3.1) becomes a kind of "free electron" model, because the kinetic term eq.(3.2), which tends 
to make electrons apart, balances with the attractive term eq.(3.3). The GWF is the very state to 
express such a situation. Thus, away from the supersymmetric case, the GWF has to be modified to 
include intersite correlations repulsive or attractive. 

The above agreement seems in close connection to the fact that the ground state of the long-range 
t-J model discussed in §4 is nothing but the GWF. 

According to the exact result, the Tomonaga-Luttinger liquid (TLL) is realized in the wide range 
of n and J/t. For the ground state, correlation functions show power-law behaviors. Hellberg and 
Mele pa] introduced a trial function ^tll which has an essentially long-range correlation factor, 



7](r) 



\N . 




— sin 


->•)! 




\N J. 



(11) 



and <3? = <£f in eq.(2.6). This factor becomes repulsive or attractive, according as the parameter 
v is positive or negative, respectively. Especially for v < —0.5 ^tll represents a separate phase. 
Correlation functions of this function show power-law behaviors, and v is analytically related to the 
TLL exponent K p by K p = 1/(2^ + 1). |26| Using this relation, we construct a phase diagram shown in 
Fig. 2, which should be compared with the diagonalization result. [27] ^tll is successful in the region 
of relatively large J/t. 

On the other hand, Fermi- liquid- type wave function ^fl 



/ s 2 r , 

r\[r) = — arctan — (repulsive 

7T C 



|22| , 23 1 with the correlation factors: 

(12) 



Fig. 1 



Fig. 2 
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en 

r)(r) = 1 H — a (attractive), (13) 

and <5 = does not show the power-law behaviors. For instance, there is a finite discontinuity 
in n(k). However bulk properties of ^fl quantitatively agree with those of vI/tll- The essential 
distinction between the correlation factors of ^fl and ^tll is in the long-distance part; the value of 
77(1)777(00) is finite for every parameter set of ^FL) but not for the other. This fact and a numerical 



experiment, [23] in which we use hybrid r/(r)'s of ^tll and ^fl connected at a certain value of 
r, indicate that the long-distance part of the correlation factor determines the long-range behavior, 
while the short-distance part does the bulk properties like energy and the amplitude of the correlation 
functions. 

For the regime of J/t ~ 0, a wave function of the form ^ = x^SF is desired on the analogy of the 
exact eigenfunction for J/t — > 0. [pg] Here x is the spin part of the wave function and i^sf is the spinless 
fermion state. Such a spin-charge separation is lacking in ^tll and \£fl- Actually, quantitatively 
reasonable results have been obtained by the form x^SF ■ fH> |29| 

In the region of small n and J/t > 2, a spin-gap state is expected. [|?]] We have not confirmed a 
corresponding state in ^tll and ^fl- On the other hand, Chen and Lee introduced a trial state for 
a gas of singlet pairs, and showed that there is a region where this function is stabler than ^tll- |^0| 
Later, the spin-gap state has been found by using Green's function Monte Carlo method. [[3l]] Fig. 3 

Finally we discuss susceptibilities of charge Xc and spin Xs, in connection with the Mott transition 
near the half-filling. |23| In Fig. 3 we show Xc calculated with ^tll- In high density region, Xc for 
every value of J/t is divergent as n — ► 1. For J/t < 2 this divergence is fitted as Xc °c 1/(1 — n). This 
divergence is due to the strong correlation effect, in comparison with the non- interacting case. In the 
meantime, Xs converges to a finite value as n — ► 1, as long as J/t > 0. These behaviors are consistent 



with the exact results for J/t = and 2. [32] In comparison with these results, the Brinkman-Rice 



transion [K3] shows different behaviors, namely Xc remains finite and Xs diverges as U — > U c (n = 1) 



or n 1 (00 > U > U c ). [34] 



3.2 Case for square lattice 

Although a lot of studies have done in 2D in connection with the high-T c superconductivity, there 
are not so many definite results. It is likely in 2D that ordered phases exist and the phase separation 



spreads to the region of relatively small J/t and n ~ 1, |35|, g6[ etc. 

First we discuss briefly the half- filled case (the Heisenberg antiferromagnet), where the ground 
state has an antiferromagnetic (AF) order. |J7| It was also shown that there exist non-AF states 



which have very close energy to the ground state by using RVB-type wave functions, [38, 39] However, 

this type of functions have not yet applied to the case n < 1. Fig. 4 

Now, let us look at the case of n ~ 1 and J/t ~ 0. In the early stage of VMC studies, the stability 
of the GWF and various ordered states is discussed. |4^, ^] They unanimously concluded a 

superconducting state with d a .2_j / 2-wave symmetry is the most stable in this region. First, we consider 
the GWF. In Fig. 4 each component of energy is shown. |4(| A conspicuous aspect different from the 
ID case is the behavior of Ej/J as n — > 1; Ej/J oc (1 — n) 0,7 , in contrast to Et/t, E3/J3, which 
behave linearly. This fact means that the GWF is unstable in itself against the phase separation for 
J/t > 0, because d 2 E/dn 2 < 0. This instability is intrinsic in the 2D GWF, and cannot be observed 
in the stabler states, as will see below. 

Next, we discuss ordered states, namely an AF state ^af: v(. r ) = 1, ^ = ^af, and a supercon- 



ducting state ^sc : r /( r ) = 1, $ = ^bcs- Here $af is a Hartree-Fock-type AF wave function, [40] and 
^bcs is the BCS-type function, for which we take different gap symmetries. In the following, *$>sc 



is used by fixing the electron number, [41, 42] which agrees with the method of the grand canonical 
ensemble [43] in the thermodynamic limit. Fig. 5 

Figure 5 shows the energy expectation values of these states for J/t = 0.5 as a function of n. As 
mentioned above, the GWF is unstable near n = 1, while ^af and ^sc of the d x 2_ y 2 symmetry are 
not. Except for the vicinity of the half-filling, d x 2_ J/ 2-wave superconducting state is the most stable. 
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For s-type and d xy symmetries, no energy reduction from the GWF is obtained in this parameter 
regime. These calculations are consistent with the previous works. 4C, 41, 43] 



Other symmetries of ^sc may be possible. Actually, we have found that many mixed-symmetry 
states including s+id are degenerate with the d rr 2_ y 2-wave state in the half- filling. Detailed studies 



on this problem are needed for n < 1. [44] 



The ^tll is extended to 2D by Valenti and Gros. |4§] They found ^tll shows power-law behaviors 
for correlation functions also in 2D. However, the energy lowering and the critical exponent obtained 
are very small, as compared with the ID system. Furthermore, since the variational energy of this type 
of functions is primarily determined by short-distance correlation factor, we can obtain comparable 
stability by ^fl- Thus a different approach may be needed to confirm the realization of the TLL state 
in 2D. g§ 

For low electron density, accurate properties are relatively easy to obtain. For the super symmetric 
case, as in ID, the GWF is stable against both ^sc and ^tll approximately for n < 0.16. However, 
in 2D this "free electron" state is unstable against ^sc with d x 2_ y 2 symmetry for higher density. 

Recently Hellberg and Manousakis have obtained a phase diagram in low electron density by 



solving few-body problems. 47] They found that in the low density limit there exists the region 



2.0 < J/t < 3.4367, where a singlet-pairing state is the most stable. Prior to this study Dagotto et 



al. found ^sc with s-wave symmetry lowers the energy around this region, [48] which is not the case 
with the ID system. At any rate, the energy lowering by \&sc an d ^tll is so small that it is not 
easy to determine the definite phase boundary in this region, except for the boundary to the phase 
separation. 

Lastly, we mention the effect of 7i^. Broadly speaking, this term contributes favorably to the 
s-type wave in the low density, but unfavorably to the d x 2_ y 2 wave for every density, when we take 
the GWF as the standard. Hence, a ferromagnetic (J3 < 0) coupling stabilizes the d x 2_ y 2-wave 
sup er conductivity. 

Details of our 2D studies will be published elsewhere. [^] 

4 One-Dimensional Supersymmetric t-J Model with Long-Range 
Coupling of Inverse Square 

It was shown that the model of the section title is exactly soluble and the ground state is the GWF. [50 
The model is written as 

ft = ZJ H(4a C J> + H - C + J ij ( S i • S i - \ n i n j)] > (14) 

i<j a 

where we require ty = Jij/2 = t/[(N/ir) sm{n x (x, — Xj)/N}] 2 }, with t > 0. It has been discussed that 
the model has a hidden symmetry called the Yangian [|5l| which explains the unexpected degeneracy 
in the spectrum. Nevertheless, explicit description of the wave functions has not been completed. 

One can use the VMC technique in judging whether a trial function is an exact eigenfunction 
or not by checking the variances of total energy in VMC sweeps, as described in §2. In this way 
we have constructed some low-lying excited states in terms of the Gutzwiller-type projection of the 



noninteracting states of up- and down-spin electrons. [53] Before we present new results, let us explain 
the notation and some previous results. 

A projected determinantal wave function V^({k}), which is our trial state, is specified by the 
occupied /c-points in the free electron state $({k}) from among as many configurations as atCV t ~nCn 1 , 
where N a denotes the number of electrons with spin a. In calculations in this section, we use systems 
of N = 41 + 2 (I: integer) and N a being odd integer with the periodic boundary condition. For 
convenience, we divide these functions into two classes: 1) states with continuous occupation of k- 
points and 2) otherwise. The class 1) has been studied, J53| but the class 2) remains to be investigated. 
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The states in the class 1) can be described using N a and current J a which is taken to be even 
integer. For the state with current J a , we construct <3>({k}) in which electrons occupy the states with 
— &f<t + irJ a /N < k < kp a + TrJ a /N. Here the Fermi wave number /cpo- is given by kp a = ir(N a — 1)/N. 
The resultant free-electron state is represented by <3?f (A^ , Jy Ny Jj_) . Thereby a generalized Gutzwiller 
wave function is obtained as ^q(N^, ^T'-^l' ^i) = 'P^F^y Jy, Ny Jj_). On the other hand, the class 
2) consists of a variety of states, which are generally not easy to be specified concisely. However, as 
will be seen below, the eigenstates in this class are restricted to the states in which only one electron 
or two are excited from the states in the class 1). 

Let us first review the case in the class 1). In this class, basic states are those without currents; J-f = 
J I = 0. In this case vI/q is an eigenfunction for every allowed values of the density n = (iV-f + N^)/N 
and the magnetization m = (jVj — N\)/N. By investigating the dependence of the energy E]\f(m,n) 
on N, m and n, it is found that the numerical results are fitted completely well by the "experimental" 
formula: 

ggftglg) = - — (Iml 3 - 3m 2 + n 3 - 3n 2 + 4n) + -^{\m\ +2n- 3), (15) 
t 12 3N 2 

which has now also been obtained analytically. |)|J] The main term was first obtained by the asymptotic 

Bethe Ansatz method, [gg] Fig. 6 

Next we consider results for finite currents. We find that the range of ^g being an eigenfunction 

is restricted by certain condition for N a and J a . In the following, we choose a typical system: n = 0.5 

and N = 60. In Figs. 6(a) and (b) the circles (filled or empty) show that ^G(Ny Jy Ny Jj) is an 

eigenstate at these values of N a and J a - The states without circles are found not to be eigenstates. 

From these results we obtain the "experimental" condition for ^g(-^t» ^T'-^i' ^i) ^° ^ e an eigenstate: 

\Ja\<N a + l, (16) 

together with 



< i—i + (17) 

2 - 2 y J 

The critical value of current given by eq.(4.3) corresponds to a momentum distribution with spin a 
where the bottom of the band becomes empty by the momentum shift. On the other hand eq.(4.4) 
suggests that a similar band exists for fictitious particles which are responsible for spin excitations. 
An exception to the condition appears in the special case of N a = 1. (53|] 

By calculating the energy increment AE caused by spin and charge currents, we find that it is 
fitted excellently by the formula 



AE / 7T 



- 



) 2 [(l-n)J c 2 + (l-m)J s 2 ], (18) 



tN V2iW 

where spin and charge currents are given by J s = (Jj — Jj)/\/2 and J c = (Jj + Jj_)/\/2, respectively. 



The formula applies regardless of values of N a , and now has also been obtained analytically. [54] 

Now we proceed to the states in the class 2). To search eigenstates from a huge number of candi- 
dates in this class, we first check small systems (N e = 6,10) thoroughly, then confirm the regularity 
for larger systems (N e = 30, etc.). As a result, we find the following facts. Eigenstates are necessarily 
the states in which at most two electrons are excited from one of the states in the class 1). Possible 
excited states as eigenstates are dependent on the values of N a and J a in the original state in the class 
1). 



On account of the limited space, here we mention the cases except for N a = 1 [53] and n = 1. [56 
In fact, the excluded cases have more abundant eigenstates probably due to their higher symmetry; 
we will report them elsewhere. In the present cases, allowed states are further restricted to those with 
only one excited electron. 

Before going to the results, let us see the representation of excitations appearing in the discussion 
below. In these excitations, one electron situated at some /c-point in the continuous occupation of the 
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original class- 1) state is excited to the nearest unoccupied Appoint. In particular, the electron at the 
middle of the cluster may move to both sides. This type of excitation is indicated by I and spin a as 
in Fig. 7. Fig. 7 

First we consider the non-magnetic case. In Fig. 6(a) possible excitation from each original class- 1) 
eigenstate (mother state) is summarized for the same system. The mother states with empty circles 
have no excited eigenstate (daughter state) from them, and ones with filled circles have at least one 
daughter state. L-shaped boxes with numbers of I show that mother states in them have daughter 
states of the specified I. A state with an arrow means only the spin specified by the arrow can excite, 
and for a state without an arrow, electrons of either spin direction can excite. Therefore, for instance, 
for the mother state (Jf, Jj_) = (—10,-12) there exists a daughter state indicated by [2,j]. For the 
states with i = —8, —7, 7 and 8, the situations are a little complex; for the mother state (0,2) there are 
two daughter states [8,f] and [—7, J,]. Moreover, four daughter states belong to the ground state (0,0). 
However these excitations are high energy processes, so that low-energy particle-hole excitations from 
the ground state cannot be described by this type of wave functions. 

Next we take up the cases with finite magnetization. In Fig. 6(b) possible excitations are depicted 
similarly. In this case, since available spins for excitation are common in the rectangular box of £, we 
show arrows together with I. The low-energy excitations are prohibited from the ground state (0,0) 
also in this case. 

To summarize, we have applied the VMC method to identify considerable amount of exact excited 
states, which are described by generalized Gutzwiller functions. These calculations are a kind of 
"experiment" . Part of the eigenfunctions are consistent with the recent analytic theory. |57| However, 
some eigenfunctions in the case of n = 1 are unexpected. We hope that the VMC results will stimulate 
further analytic theory on this supersymmetric model. 



5 Summary 

In this paper we have described recent developments and new results for t- J-type models by the VMC 
method. As we have seen above, as well as in refs. 7 and 8, the VMC method have contributed to 
various issues of the strong correlation. Keeping its merits and demerits in mind, one can develop it 
further also in other systems of interest, for example, the multiband models [58, 59] and the fractional 



quantum Hall effect. [60] 
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Fig. 1 

Comparison of the total energy per site between the GWF and Bethe ansatz [24] and the non- 
interacting system for J/t = 2 as a function of electron density. Although the VMC results for 
N = 102 is plotted, the analytic expression is available. |23|] 
Fig. 2 

Phase diagram of the ID t-J model obtained by the ^tll- The curves show the contours of constant 
correlation exponent K p . The used system has 100 sites. 
Fig. 3 

Charge susceptibility in unit of 1/t vs. n for some values of J/t. Symbols are the results of the 
vE'tll- Solid lines for J/t = and 2 |32[ ] represent the exact analytic values. Dashed line is for the 
non-interacting system. Dotted lines for J/t = 1.0 and 2.5 are a guide to the eyes. The sizes of the 
symbols represent the relative amplitude of possible error. 50 ~ 210-site systems are used. 
Fig. 4 

Expectation values of each energy components E t = {Tit), Ej = (Wj), E3 = as a function 

of electron density. The used systems are TV = 10 x 10 ~ 26 x 26 with the periodic (x axis) and 
antiperiodic (y axis) boundary condition. 
Fig. 5 

Comparison among some variational energies as a function of electron density. The used systems are 
N = 10 x 10 (diamond), 12 x 12 (square), 14 x 14 (circle), 16 x 16 (downward triangle), 20 x 20 
(upward triangle) with the same boundary condition with Fig. 4. 
Fig. 6 

Condition of N a and J a for being an eigenstate in the class 1). The states specified by solid and open 
circles are eigenstates. Two typical cases of different values of m are shown: (a) non-magnetic case, 
(b) finite magnetization. The size of the system and the electron density are fixed at N = 60 and 
n = 0.5, respectively. The difference of solid and open circle, boxes, arrows and numbers indicate 
possible excited eigenstates in the class 2) constructed from the original eigenstate in the class 1); see 
text and Fig. 7 for detailed explanation. 
Fig. 7 

Representation of excited states (class 2)) from an original state in the class 1). Shown is only the 
o"-spin configuration, which varies in the excitation, in fc-space. The configuration of — a-spin does 
not change. An solid (open) circle represents an occupied (unoccupied) Appoint. [£, a] denotes the 
following excitation; the \£\-th (1 < \£\ < (N a + l)/2) electron (spin a) from the nearer boundary, 
namely k = sgn(£)/cFo- + ^J^/N, is displaced to the nearest unoccupied fc-point from the original 
class-1) state, which is specified by J a and N a , as in Figs. 6(a) and (b). 
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